Metallic Nanosphere in a Magnetic Field: an Exact Solution 
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We consider the electron gas moving on the surface of a sphere in a uniform magnetic field. 
An exact solution of the problem is found in terms of oblate spheroidal functions, depending on 
the parameter p — "I>/$o, the number of flux quanta piercing the sphere. The regimes of weak and 
strong fields are discussed, the Green's functions are found for both limiting cases in the closed form. 
In the weak fields the magnetic susceptibility reveals a set of jumps at half-integer p. The strong 
field regime is characterized by the formation of Landau levels and localization of the electron states 
near the poles of the sphere defined by a direction of the field. The effects of coherence within the 
sphere are lost when its radius exceeds the mean-free path. 
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The electronic propecl;ies of cylindrical and spherical 
carbon macromoleculesEJ have attracted much of theo- 
retical interest last years. A large part of these stud- 
ies is devoted to the band structure calculationsB and to 
the effects of topology for the transport and mechanical 
properties of these nanostructures. The main interest 
to the topological aspects is connected with the carbon 
nanotubes, where one can investigate the effective mod- 
els, based on the band calculations^nd incorporating the 
particular geometry of the object. tl 

At the same time, the problem of the topology appears 
to be a general one and may be studied independently 
from the physics of fiarbon materials. The recent ad- 
vances in technologyQ let one think of a wider class of 
the sphericsd nanostructures, e.g. the spheres coated by 
metal films,u whose properties may differ from those of 
planar objects. 

In this paper we study the gas of electrons moving 
within a thin spherical layer in an applied magnetic field. 
We find the exact solution of this problem in terms of 
oblate spheroidal functions. This physical application of 
the theory of spheroidal tinctions is demonstrated ap- 
parently for the first timcBia We show the jumps in the 
susceptibility of the system in "weak" fields and the local- 
ization of the electronic states in "strong" fields. The last 
effect could be experimentally investigated for the hemi- 
spherical tips of nanotubesE At the intermediate fields 
the sophisticated structure of the functions makes the an- 
alytical treatment impossible and the numerical methods 
should be used for the analysis of observable quantities. 



I. A SETTING-UP OF THE PROBLEM 



Let us consider the electrons moving on a surface of the 
sphere of radius rg ■ In the presence of the uniform mag- 
netic field B we choose the gauge of the vector potential 
as a vector product A = i(B x r). The Hamiltonian of 
the system is given by 



n 



2m, 



■(-^V + eA)2 + [/(r), 



(1) 



where me is the (effective) mass of an electron, — e its 
charge; we have set h — c = I and omitted the trivial 
term connected with the spin of the particle. We assume 
that the total number of particles N (with one projec- 
tion of spin) is fixed and defines the value of the chemi- 
cal potential /i and the areal density v = N / (47rrQ). The 
confining potential U{r) — Q aX tq < r < + 5r and 
U{r) —* OO otherwise. We focus our attention on the 
limit (5r ^ ro , when the variables of the problem are sep- 
arated. The radial component R{r) of the wave function 
is a solution of the Shrodinger equation with the quantum 
well potential. Henceforth we ignore the radial compo- 
nent and put r = To in the remaining angular part of the 
Hamiltonian Tiji ; it can be done if /i lies below the first 
excited level of R{r), which in turn means Sr < i/^^/^. 
We choose the direction of the field B along the z— axis 
and as a north pole of the sphere {9 — 0) and look for the 
eigenfunctions to the Shrodinger equation TCn'i' = E'^/ in 
the form ^'(0, 0) = S'(0)e""'^. Defining the dimensionless 
energy e = 2mergi? and setting rj = cos 9, we write 
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drj 
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2mp 



s = o. 

(2) 



We introduced here the important dimensionless param- 
eter 

p = eBrl/2 = ^Br^/^o = rl/{2ll) = (3) 

with the magnetic fiux quantum $0 = 27r/e = 2 • 10^^^ 
T-m^, the magnetic length /, — (eB)^^^^, the cyclotron 
frequency Uc = eB/rrie- Note that for a sphere of radius 
To = 10 nm one has p = 1 at the field i? ~ 6 T. 

The equation (^ is known as the spheroidal differs 
tial equation and was extensively studied previously.! 
The solutions to it are given by the oblate (angular) 
spheroidal functions Sim {p, ">]) with the corresponding 
eigenvalues eim{p). 



1 



It is known that spheroidal functions belong to the 
simplest class of special functions which are not essen- 
tially hypergeometric ones. For the spheroidal functions 
there are no recurrence relations, generating function rep- 
resentations etc., which are characteristic for the clas- 
sical special functions. The spectrum eim{p) is found 
as the eigenvalues of the infinite matrices. For different 
sets of functions orthogonal on the interval (—1, 1), these 
matrices are reduced to the tridiagonal form and could 
be further analyzed within the chain fractiaruformalism 
or by explicit (numerical) diagonalizationB^Q Perhaps, 
the mostly known quantum-mechanical application of the 
spheroidal functions is the problem of an electron in two- 
center Coulomb potential (Hj molecule). 



completely filled and (L-l-l)th level is empty. Linearizing 
the spectrum we have ejm ~ 2L{1 — L—l/2+pm/ L). It is 
clear that at p > 1/2 the state {I = L + l,m = —L — 1) is 
energetically more favorable than the state {I — L,m — 
L), with the resulting change in the occupation of the 
levels. 

This phenomenon is accompanied by the jumps in 
the static susceptibility. Consider the free energy F — 
-^E;mlii(e^^~'^""-'^^ + !)■ Then, apart from the Pauh 
spin contribution, the magnetic (differential) susceptibil- 
ity X — dM/dB ~ —d'^F/dB'^ is given by 



X = 



2 2 

hn 



dp^ 



-nF{eim) + 
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dp 



II. THE WEAK FIELD, SUSCEPTIBILITY 

Let us first concentrate on the case of the weak field. 
In the absence of the field, the spectrum is that of a free 
rotator model and the-splutions to (J^ are the associated 
Legendre polynomialslU: 



e;m(0) = ;(/ + !), 



121 + 1 {I - \m\) 
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(4) 



According to (^ the wave functions are the spherical har- 
monics '^irn{0,<p) = r'^^Yim{Q , (j)) and are normalized on 
the surface of the sphere, Tq J \ ^\'^ s\n9d9d(j) — 1. This 
normalization facilitates the comparison of our results 
with the case of planar geometry. 

For p ^ Q one develops the perturbation theoryj-pxpund 
the initial wave functions (^ as long as p^ < 4Z.ErEl The 
energies and wave functions in this case are given by se- 
ries in p^ : 



with the Bohr magneton plb — e/2TOe, the Fermi function 
np{e) and n'p{e) its derivative. The first term here gives 
the diamagnetic contribution x"^ — — 2/3/i^merQ and 
the second term is the paramagnetic one. At low tem- 
peratures ujc ^ T 't^ ^ we let n'p{Eim) ~ — 5(e;m — m) 
and change ~ dm. Performing the integration 
we get 

X = N\^xl/2^l) [-2/3 + Y.fi-L- 1/2)2] 

where the summation is restricted by |^ — L — 1/2| < p. 
From (^ we see that x exhibits the jumps at p = 
1/2, 3/2, 5/2 . . . while x — ^ in the formal limit p 00. 
This behavior is shown in the Fig. |], together with the 
quantity AI/B. At the other fillings, i.e. at the other 
values of fi the jumps of x take place at other values of p 
and the qualitative picture remains true. The amplitude 
of the jumps is A^ times larger than the Pauli spin sus- 
ceptibility ~ Nfi"^/ fj,; this coherent effect vanishes if the 
coherence on the sphere is lost due to the finite quasipar- 
ticle lifetime (see below). 



£irn{p) = l{l + l) + 2pm + 
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To clarify the criterion p^// < l we consider the case 
when the field ceases to be small for electrons near the 
Fermi level I ~ tqV^^, mostly contributing to the phys- 
ical properties. This case corresponds to the following 
relation 



BrI 
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Hence the field cannot be treated as a perturbation, if 
the energy of the magnetic field in the volume of the 
sphere is 10^ times larger than the characteristic plasma 
frequency. For densities v ~ 10"'^^ cm"^ and tq = 10 nm 
it corresponds to the fields greater than 40 T. 

For the weak field regime it is interesting to observe the 
following property of the spectrum (^ . For simplicity we 
consider the situation when the Lth unperturbed level is 



III. THE WEAK FIELD, GREEN'S FUNCTION 

Let us discuss the properties of the electron Green's 
hmction. For the two points on the sphere r ^ [9, 0) 
and r' ^ (9' 6) we define 



G(r,r',c.) = ^ 



U + PL- Elm 



(8) 



hn 

with Eijn = (2mer^)"^e/m. 

In the absence of the magnetic field ^'/m(^, < 



^0 ^Ira ( 



and one can find (see below) an exact rep- 



resentation of G through the Legendre function 



2 cos na 



P-l/2+ai-COSn), (9) 



where we introduced a = ■y/2mer§(/x -|- oj) + 1/4 and the 
distance ft on the sphere : 
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COS r2 = COS 6 cos 6' + sin sin 0' cos 

We immediately see from the logarithmic singular- 
ity in the one-point correlator (fi = 0) as it should be. 
It is instructive to find out how the usual expressions 
for the planar geometry are recovered for the large ra- 
dius of the sphere, tq — > oo; in this case a oc rp, while 



r — r 



/ro = r/ro. 



In the limit a 3>i-4i and for 

,-1 .to 



G"(c.) 



cos(a7r — aVl — it /A) 



\/27ra sin f2 



cos 7ra 



(10) 



The existence of two oscillating exponents exp±i(a7r — 
aVL — 7r/4) corresponds to quantum coherence of two 
waves. One propagates along the shortest way between 
two points and another wave goes along the longest way, 
turning round the sphere. In the theory of metals one 
considers jwl <C /i = kp/{2mf,), while the finite quasipar- 
ticle lifetime can be modelled by ascribing the imaginary 
part to to. We write in this sense ^/2rn^ = kp + i/lmfp 
with the mean free path Imfp- When tq ^ Imfp, we see 
that the coherence breaks and the only surviving expo- 
nent in ([l0| ) has the form 



exp 



ikpr ~\- i- 



I'mfp 



(11) 



which expression coincides with the usual findings .0 

It is possible to obtain the closed form of the Green's 
function in the weak field regime. We sketch the corre- 
sponding derivation below. 

First we neglect the p^/l terms in the wave func- 
tion (|) and write YC^{e' ,^)Yi^{e,Q) = (4^)^i(2? 
l)e-""'*P;"(cos6')Pf'"(cos6l'). Using (|) and dropping 
the terms containing p^/a ~ /I, we have 



21 + 1 
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pm^ , 1 
a 2 



Now we represent the appearing fraction as the Taylor 
series ia+pm/a— z)^^ — exp (^^) (a— z)^^ and sub- 
stitute m by in the last exponent. After that we can 
sum over m m (1)0 and represent the remaining sum 
over I as the contour integral, to obtain 



G(u;) 



exp 



exp 



.p 



a dadcf) 
,p 3^ 
a dadcf) 



dv 



Pu{— cosil) 



2m 2 sin irv a ^ v ~ 1/2 
G"iu;), (12) 



with G°(w) given by (|^). When a 3> 1 the expression ( p^ ) 
could be simplified at the above condition a sin f2 ^ 1. In 
this case we use (10) and observe that, upon differenti- 
ating over 0, the main contribution of order of a stems 
from the numerator exp(±iaf2) of dlQ). Then one uses 
the identity exp [f ^a] = ^ exp [z^T^'* ^^'^ finds 



TOe exp[ipP(TT — fl)] 

2 \j2'Ka sinfi 

g — m(7r — n)+27r/4" 



pia(7r— n) — i-7r/4 



COS TT a ■ 



-pH) 



cos 7r(a — p/3) 



(13) 



with /3 = sin ^^ sin 6*' sine/)/ sin ri. As before, in the pres- 
ence of damping tq > Imfp^ the coherence on the sphere 
breaks and the only surviving term in (|l^) has the form 

Gda,nped{i^) ^ exp(iz(r' X r)/2/2)G;]„„pgrf(w), (14) 

in accordance with previous results (see e.g.0). 

IV. THE STRONG FIELD REGIME 

Let us now turn to the case of strong fields, p — > oo. We 
define the integer number n > as follows : 2n — I ~ \m\ 
for even I — \m\ and 2n-f 1 = ^ — |m| for odd I — \m\. The 
value of n has a simple meaning, it corresponds to the 
number of zeroes of wave function Sim (p, cos 9) within the 
interval 9 £ (0, 7r/2), by analogy with the known prop- 
erty of (cos 6*). The spectrum of (||) is then given by 
a series 



eimip) = 4p [n + (m + |m| + l)/2] - (s^ 
+0 [s^p] , 



m 



1)12 
(15) 



with s = 2n + |m| + 1 ; one has s = Z + 1 {s = I) for even 
(odd) values of Z — to. The eigenfunctions are given by 



^hn V) = Sim (P, V)±Slm (p, -7?) , 



(16) 



where plus (minus) sign corresponds to even (odd) values 
of I — niB The functions Sim{p,v) £^re found as a series 
in the Laguerre polynomials L^{x) ; in the main order 
of they can be written as 

Simip, v)^cil- ,7)l™l/2e-^'(i-'') 4^1(2^(1 - 77)), (17) 



2|m|p|m| + l^| 



27rr2(n -|- |to|)! 



1/2 



We see that in the main order of p the Landau quantiza- 
tion takes place, i.e. the spectrum ( p^ is that of quantum 
oscillator jjsiith the cyclotron frequency being the energy 
quantum. tj The convergence and hence the applicability 



of the series (15) is. given by the condition p ^ s ~ /. In 



its turn, it meanslij that all n zeroes of the approximate 
eigenfunction Lj™'(2p(l — ry)) in ( p^ lie in the northern 
hemisphere, ry > 0, as they should. 

We notice the following important property of the spec- 
trum (p5|). For given non-positive to the values of 
corresponding to Z = 2n -I- |to| and / = 2n -f |to| -f 1 co- 
incide. This property and the form of the wave function 
( p6| ) can be understood as follows. At p ^ oo the field- 
induced potential sin^ 9 in (0) localizes the particles 



3 



near the poles of the sphere. This form of two-well po- 
tential leads to the discussed degeneracy of energy levels, 
while the total wave function is given by a symmetriza- 
tion (16) of the wave functions related to each of 
the wells. The possibility of quantum tunneling between 
the wells lifts the degeneracy and produces the exponen- 
tially small energy splitting £r^ e~^P) between the states 

It would be tempting to obtain the Green's function 
at p — > oo in the closed form. It can be done at the 
following simplifications. First, we neglect the expo- 
nentially small splitting and notice that 5+ (77)6'+ (77') + 
S-{ri)S-{r]') = 2[S{r])S{r]') + S{-ri)S{-r)')l i.e. the cor- 
relations within one hemisphere only survive. Now we 
leave only the leading terms (|l7| ) of the wave functions 
and 0{p) terms in energies (ISp. The necessity to re- 
strict the summation in (j^) by the terms with s p 
could be modelled by inclusion of the cutoff factor e^'** 
with S p~^ into (^). We put /i = for simplicity 
of writing and raise the denominator into the exponent, 
(uj - Eim)~^ = i /o°° e^*(^""-"). Next we perform the 
summation over n with the use of bilinear generating 
function for the Laguerre polynomialai3 : 



(n + \m\)l 

n—O ^ 11/ 



4"i(2;)4r'(2/)^" 



(xyz) 



-|m|/2 



1 



■ exp — z 



'|m| 



with the modified Bessel function Imiw); in our case 
X = 2p{l ~ rj), y = 2p{l - rj') and z = e'*""^. Mak- 
ing use of the property — i™'J„i{w) and rescal- 
ing tuJc/2 — > i we arrive at the intermediate formula for 
77,77' > : 



dt 



An J sin(t -I- i5) 



m— — oo 



xy 



sm{t + id) 



The summation over m is now easily done 
(Eme™'^J™(w) = e™'''"'^) and we obtain the Green's 
function as a sum of two terms, referring to two hemi- 
spheres. The "northern" term, corresponding to the 
above expression, is given by 



i(5-t-oo 



[uj] ~ — - — e / -; — e 



(18) 



i5 



ir/2+i(5 



rriee 



Sttcostt- 



dt- 



cos t 



-Tr/2+ii5 



with the analogue of the distance on the sphere p — 
2p[2 - 7? - 7?' - 2^(1 - ?7)(1 - 7;')cos(0 -I- iS)] and that 



of the vector product v = 2p^{\ — ri){\ — 77') sin((^-|-i5). 
The form of the "southern" term is obtained by chang- 
ing T] ^ —7] and 77' —77' in these expressions. The last 
integral is reduced for vanishing S to the hypergeometric 
function "^{a, b; z) and we obtain the "northern" term in 
the form : 



G"(w) 



V— p/2 — 5uj / uJc. 



An 



00 



Lnip) 



1/2- 



(19) 



The last equation is similar to previous finding 

Let us discuss the applicability of (|l9|). Our deriva- 
tion was straightforward until (|l^), while the last step 
demanded S 1/p 0. The incomplete restructuring of 
the spectrum into the Landau level scheme is absent in 
the usual planar geometry, wherein one would put 5 = 
and the expression ( p^ ) would be exact. In the spherical 
case we cannot treat the higher levels with / > p in an 
analytical way, it is mimicked in ( |T9| ) by the appearance 
of the exponential cutoff at a; ptOc- 

A subtler issue in justifying ( p^ ) is the shift of the last 
integration in dH) to the interval (-7r/2, 7r/2). The inte- 
gration over the remaining segments (±7r/2,±7r/2 -I- iS) 
yields the factor cos(7ru;/a;c), thus the contribution of 
these segments to the Green's function has no poles in 
Lu. This contribution could be combined with the smooth 
(real) part of G(a;) stemming from the consideration of 
the highest energy levels. 

As a result, one may conclude that the expression ( pj| ) 
correctly reproduces the basic properties of the Green's 
function for uj < puJc- 

The finite value of the cutoff parameter S in ( pj| ) be- 
comes important for the one-point correlation function, 
i.e. at rj = rj' and cf) = 0. The residues of the Green's 
function define the local density of states (LDOS) by the 
relation 7V(r) = / dcj np{uj)[G{r,r,uj — iO) — G{r,r,uj + 
iO)]/ (27ri). If we assume that n lowest Landau levels are 
filled by the electrons, then it follows from (^9|) that the 
LDOS is given by N{r) = npe'^'-P/^ /{2nr'^). In this case 
the finite S ^ 1/p provides a smooth variation of LDOS 
of the form 

N{r) cx exp[-0(sin2 6l)], 

which variation is absent in the usual planar geometry. 
This result, however, may depend on the approximations 
made and requires further numerical investigation. 

It should be stressed that at the intermediate fields, 
at p > 1 and p ^ I ^ p^, the spectrum and wave func- 
tions are not principally reduced to closed expressions 
of hypergeometric type. The numerical methods are in- 
dispensable here. We calculated the evolution with p of 
the energy levels with I = 0, ... 5 by diagonalizing the 
200 X 200 tridiagonal matrices in the basis of P™{ri). 
The results are shown in the Fig. ^ where we made a fol- 
lowing convention. From the general property of Eq.(0) 



4 



it foUaajs that ei^-m{p) = £i,m{—p) with formally nega- 
tive p.O For better readability of the graph, it is possible 
then to show all the data, plotting only half of them with 
one sign of m but in the formally extended region of p. 
One can see in the Fig. || that the degeneracy at p — 
is eventually changed by the Landau levels formation at 
p = 10. In the intermediate region p ~ 3 the absence of 
any structure in the levels' scheme is noted. The mesh 
of lines at the intermediate fields mimicks the chaotic 
behavior, although this is not so. Both in the quantum 
problem considered here and in its classical coutcrpart 
one has two variables (9, (f) and two integrals of motion, 
the energy and the projection of the angular momentum 
onto the field direction. To obtain chaos, it suffices to 
break the rotation symmetry, — > + 5(j). The latter 
problem is, however, beyond the scope of this study. 

In conclusion, we demonstrate the exact solution of the 
electron gas on the sphere in the magnetic field. In the 
limits of weak and strong fields this solution is reduced to 
the hypergeometric functions and the observable quanti- 
ties are found in the closed form. In the case of intermedi- 
ate fields, the solution is not essentially hypergeometric 
function and the observables require further numerical 
analysis. 
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FIG. 1. The dependence of the static susceptibility x on 
the number p of magnetic flux quanta piercing the sphere. 
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FIG. 2. The dependence of the energy levels eim on 
the number p of magnetic flux quanta piercing the sphere. 
For convenience of presentation, we plotted the evolution of 
Eim{p) with m < on the rhs of the plot, at p > 0. The 
evolution of eimip) with m > is depicted on the Ihs of this 
plot for the formally negative p, see the text for additional 
explanations. 
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